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ABSTRACT 

The fine-structure lines of singly ([Ne II] 12.8 /im) and doubly ([Ne III] 15.6 /mi) ionized neon are among the 
most prominent features in the mid-infrared spectra of star-forming regions, and have the potential to be a powerful 
new indicator of the star formation rate in galaxies. Using a sample of star-forming galaxies with measurements of 
the fine-structure lines available from the literature, we show that the sum of the [Ne II] and [Ne III] luminosities 
obeys a tight, linear correlation with the total infrared luminosity, over 5 orders of magnitude in luminosity. We 
discuss the formation of the lines and their relation with the Lyman continuum luminosity. A simple calibration 
between star formation rate and the [Ne II] +[Ne III] luminosity is presented. 

Subject headings: galaxies: ISM — galaxies: starburst — infrared: galaxies 



1. INTRODUCTION 

The star formation rate (SFR) is one of the most fundamen- 
tal parameters in galaxy evolution, but not trivial to measure. 
Practical methods of estimating the SFR ultimately rely on in- 
direct indicators of the hydrogen Lyman continuum luminosity. 
From this, one estimates the number of massive, ionizing stars, 
and for an assumed initial mass function, the total stellar mass. 
The most commonly used SFR estimators make use of the op- 
tical hydrogen recombination lines and [O II] A3727 (see, e.g., 
Kennicutt 1998 for a review). Depending on the application, a 
variety of estimators based on continuum emission have been 
devised and intercalibrated, including X-rays, mid-ultraviolet, 
far-infrared (IR), and radio synchrotron (e.g., Bell 2003; Gil- 
fanov et al. 2004). 

The advent of the Infrared Space Observatory (ISO) and es- 
pecially the Spitzer Space Telescope has brought renewed inter- 
est in mid-IR diagnostics of star formation. Two mid-IR SFR 
estimators have been discussed: the strength of polycyclic aro- 
matic hydrocarbon (PAH) features (Peeters et al. 2004; Wu 
et al. 2005) and the mid-IR continuum (Roussel et al. 2001; 
Wu et al. 2005). These estimators, however, are not without 
complications. The PAH features are known to be exception- 
ally weak in low-metallicity galaxies (e.g., Wu et al. 2006, 
and references therein), their strength being governed by com- 
plicated, but poorly understood formation and destruction pro- 
cesses (O'Halloran et al. 2006; Wu et al. 2006). PAH emis- 
sion arises from a variety of interstellar environments, not only 
star-forming sites, and thus its ability to reliably trace starburst 
regions remains a matter of debate (Haas et al. 2002; Peeters et 
al. 2004). The mid-IR continuum strength, on the other hand, 
shows considerable variation in its relation to the total IR lu- 
minosity for star-forming galaxies (Dale et al. 2005; but see 
Takeuchi et al. 2005). Moreover, without additional diagnos- 
tics it is often difficult to tell whether the primary source of 
dust heating responsible for the mid-IR emission derives from 
star formation or active galactic nuclei. 

This paper proposes a new SFR indicator based on the 



strength of the fine-structure transitions of [Ne II] A12.814 /im 
( 2 P\/2 P3/2) and [Ne III] A15.554 /im (M -^ 3 P 2 ). [Ne II] 
is an excellent tracer of ionizing stars because (1) its ioniza- 
tion potential (21.56 eV) is higher than that of hydrogen; (2) 
it is a dominant ionization species in H II regions, being one 
of the principal coolants (Burbidge et al. 1963; Gould 1963); 
(3) its relatively high critical density (n c = 4.3 x 10 5 cm" 3 ; Pet- 
rosian 1970) ensures that the line flux is insensitive to density, 
since most H II regions have electron densities substantially be- 
low n c ; (4) it is abundant (solar Ne/H = 1.20 x 10~ 4 ; Grevesse 
& Sauval 1998), with a nucleosynthetic history closely follow- 
ing that of oxygen (Garnett 2002); and (5) its long wavelength 
significantly reduces its sensitivity to dust extinction compared 
to optical lines (the extinction at mid-IR wavelengths is only 
a few percent of the extinction at optical wavelengths). Dou- 
bly ionized Ne shares these properties, with the difference that 
[Ne III] can be the dominant ionic species in regions of lower 
density (n e < 100 cm" 3 ) and higher excitation (blackbody ra- 
diation temperatures T r > 4 x 10 4 K). For example, low-mass, 
low-metallicity galaxies, which are known to exhibit systemat- 
ically higher excitation, show a greater fraction of ionized neon 
emission in the [Ne III] 15.6 /im line rather than in the [Ne II] 
12.8 /zm line (O'Halloran et al. 2006; Wu et al. 2006). 

We demonstrate, using a sample of star-forming galaxies 
with mid-IR spectra, that the sum of [Ne II] + [Ne III] lumi- 
nosity tightly correlates with a well-known measure of star for- 
mation, the total IR luminosity. We discuss the ionization of 
Ne and the connection between the [Ne II] and [Ne III] emis- 
sion and Lyman continuum luminosity, from which a simple 
relationship between [Ne II] + [Ne III] luminosity and SFR fol- 
lows. 

2. DATA 

To explore the utility of the mid-IR fine-structure lines of 
Ne as a SFR indicator, we have assembled from the literature 
data for star-forming galaxies with measured [Ne II] and [Ne III] 
fluxes. We have strived to ensure that the data were obtained 
with a sufficiently large aperture so that most of the emission is 
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included. Since the [Ne II] line sits adjacent to the prominent 
12.7 /j,m PAH feature, it is also important that the spectra have 
sufficient spectral resolution so that the line can be unambigu- 
ously resolved. Lastly, while the resulting compilation is neces- 
sarily heterogeneous, we have tried to minimize systematic ef- 
fects by concentrating only on a few relatively large, uniformly 
analyzed samples from the literature that were observed with 
the same instrumental setup. The final sample of 57 sources, 
summarized in Figure 1, contains 34 galaxies observed with the 
Short Wavelength Spectrometer (SWS) on ISO (Genzel et al. 
1998; Thornley et al. 2000; Verma et al. 2003) and 23 galax- 
ies observed with the Infrared Spectrograph (IRS) on Spitzer 
(O'Halloran et al. 2006; Wu et al. 2006). Although some of the 
sources are well-known composite systems that contain both 
an active galactic nucleus and a starburst (e.g., Arp 220 and 
NGC 6240), it has been previously established that the bulk of 
the bolometric luminosity in these sources is dominated by the 
starburst component. Our diverse sample spans a wide range of 
galaxy types, from low-mass, low-metallicity dwarf galaxies to 
massive merger remnants. Among the 57 objects in the sample, 
16 are early-type spirals (Sa-Sbc; (M B ) = -19.4 mag), 7 are 
late-type spirals (Sc-Sm; (M B ) = -18.8 mag), 14 are irregulars 
and blue compact dwarfs (Im, 10, BCD; (Mb) = -17.1 mag), 
and 20 are well-known merger or strongly interacting systems 
with peculiar morphologies ((M B ) = -20.4 mag). 

Figure 1 illustrates that the [Ne II] emission strongly corre- 
lates with the IR emission, over 5 orders of magnitude in lumi- 
nosity. An ordinary least-squares regression of L iR on L[ Ne h] 
yields 

legion] = (1.01 ± 0.054) logL IR - (3.44 ±0.56), (1) 

with a scatter of 0.51 dex. The correlation appears to steepen 
toward low luminosities (Li R < 10 9 5 L Q ), in that less [Ne II] 
emission is produced for a given amount of IR luminosity. 
The lower-luminosity galaxies have a higher ratio of [Ne III] 
to [Ne II] emission, suggesting that in these galaxies a sig- 
nificant fraction of the Ne is doubly ionized, most likely be- 
cause in a low-metallicity environment the ionizing radiation 
field hardens as a result of reduced line blanketing and block- 
ing in stellar atmospheres (e.g., Thornley et al. 2000; Mad- 
den et al. 2002). We examine this effect by plotting the sum 
of [Ne II] and [Ne III] (solid points in Fig. 1), for those ob- 
jects for which both lines are measured. Indeed, the lower- 
luminosity points now exhibit significantly less scatter and fol- 
low more closely a smooth continuation of the trend defined 
by the higher-luminosity systems. An ordinary least-squares 
regression of Li R on L [Ne n ]+[Ne m] yields 

logL [Ne n]+[Ne m] = (0.98 ± 0.069) logL IR - (2.78 ±0.70) (2) 

and a scatter of 0.49 dex. 

The IR luminosity is a robust indicator of the SFR in star- 
forming galaxies because a significant fraction of the bolomet- 
ric luminosity of young stars is absorbed and reradiated by dust 
grains. We expect the IR and Ne luminosities to be related be- 
cause, as discussed in the next section, the lines of ionized Ne 
directly trace the Lyman continuum, which is itself produced 
primarily by young stars. The tight correlation between the lu- 
minosity of ionized Ne and the IR luminosity indicates that the 
luminosity of the Ne fine-structure lines can be used as an ef- 
fective substitute to estimate the SFR. The theoretical basis for 
this expectation is examined in the next section. 



3. THE RELATION BETWEEN THE NEON LINES AND SFR 

3.1. The Ionization ofNe 

The relative abundances of the ionic species are set by a bal- 
ance of the rates of photoionization P; and recombination R\ for 
each species 

N i P i =N M R i n e . (3) 

The photoionization rate is the integral over a frequency- 
dependent absorption cross section a\{v) and the ionizing ra- 
diation J(i>) above the ionization threshold v{. 

Pi = / — J{v) aM dv. (4) 

In the case of Ne, Figure 2 shows that the absorption cross sec- 
tion of Ne + peaks at a frequency corresponding to a Wien tem- 
perature of 1.5 x 10 5 K, which is higher than the characteristic 
radiation temperature of early-type stars (T r w 4 x 10 4 K). Thus, 
in a region ionized by typical massive stars the dominant ionic 
species would be Ne + . However, as shown in Figure 3, [Ne III] 
will begin to replace [Ne II] in ionized gas that is at lower den- 
sity or subject to intense or high-temperature radiation. The 
ionization equilibrium shown in Figure 3 is calculated in the 
absence of a specific model by expressing the mean radiation 
as diluted blackbody radiation, 

AitJ v = ttTB v , (5) 

where T is analogous to the r~ 2 geometrical dilution in a spher- 
ical model. In Equation 3, the dilution factor may then be com- 
bined as a ratio, T/n e , with the density, which is also unspec- 
ified in the absence of a specific model. Figure 3 illustrates 
how the use of the combination of [Ne II] + [Ne III] improves 
the robustness of Ne as an indicator of the SFR, by extend- 
ing its applicability to galaxies with higher T/n e , for example, 
to lower-luminosity, lower-metallicity galaxies that have higher 
excitation (e.g., Thornley et al. 2000; Madden et al. 2006). 

3.2. The Line Intensity 
The general expression for the intensity of a spectral line is, 

-i _? -i hv A21 f 

I (ergs s cm z sr ) = — / n 2 d£, (6) 

4tt J 

where A21 is the Einstein coefficient for spontaneous emission, 
«2 is the number density of ions in the upper state, and the in- 
tegral is over pathlength dl. The population of the upper state 
is in general found from the equilibrium set by the rate equa- 
tions for collisional and radiative transitions between all levels 
of a species. The radiative rates are considerably faster than the 
collisional rates at electron densities below the critical density 
for collisional de-excitation, ~ 10 5 cm -3 , for both the [Ne II] 
and [Ne III] fine-structure transitions. Therefore, for lower den- 
sities, it is possible to approximate the rate of emission, A 2 \n 2 , 
by the upward collisional rate alone. In the case of the [Ne II] 
doublet, the line intensity is then approximately 

I lNeU] = ^f + C 12 EM, (7) 

47T 

where 7 is the interstellar abundance of neon relative to hydro- 
gen (9.9 x 10" 5 ; Simpson et al. 1998), /+ is the fraction of Ne 
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that is singly ionized (Fig. 3), and C12 « e is the upward col- 
lisional rate (s -1 ) between the 2 Pt,/2 and the 2 P\/2 levels, and 
EM = J n 2 dl is the emission measure. 

In the case of the [Ne III] triplet, one needs to consider the up- 
ward transitions to both higher states because a collisional tran- 
sition to the uppermost 3 Po level will result in a radiative tran- 
sition to the intermediate 3 Pi level, followed immediately by a 
transition at A 15. 6 /im to the ground state 3 P2- The rate of the 
electric quadrupole transition 3 Pq — > 3 P2 is negligible compared 
to the rates of the two magnetic dipole transitions (A/ = 1). 
Thus, the intensity of [Ne III] is approximately 



hv 

/[Ne IH] = 1-1 /+2 (C12 + C13) EM, 



(8) 



where the subscripts 1, 2, and 3 on the collision rates refer to 
the levels 3 P2, 3 P\, and 3 Po, respectively, in order of increasing 
energy. 

The collision rates are related to the collision strengths $721, 
assuming a thermal velocity for the electrons, 



Ci2 = 8.63 x 10" 8 — 



10 4 K 



1/2 



(9) 



where uj\ is the statistical weight of the lower level. For the 
[Ne II] A12.8 /xm line, Q 2 i/wi = 0.0785 (Griffin et al. 2001). 
For the [Ne III] A15.6 fim line, 2 i/wi = 0.1548 and O31/W1 = 
0.0042 (Butler & Zeippen 1994). 

Integrating over the source area A = J dil and assuming 
T e 10 4 K, the sum of the luminosities of the two lines can 
be expressed as 



Magellanic Clouds come from Giveon et al. (2002), while 
the M33 sources come from Willner & Nelson-Patel (2002fl 
Three points deserve notice. First, for the Galactic sources 
both [Ne II] and the sum of [Ne II] + [Ne III] trace Bra lin- 
early. Second, for the extragalactic sources, all of which have 
subsolar abundances (especially the Small Magellanic Cloud), 
[Ne II] systematically deviates from a linear correlation with 
Bra whereas the sum of [Ne II] + [Ne III] remains linear. Fi- 
nally, the relation given in Equation 12, with /+ = 0.75 (which 
is a reasonable estimate from the peak of the curves shown in 
Fig. 3) and fi on = 0.6, extremely well describes the [Ne II] data 
for the Galactic sources. The same relation, with a small ad- 
justment in zeropoint to account for the ionization fraction, very 
well matches the sum of [Ne II] + [Ne III] for all the sources; the 
scatter around the linear relation is only ~ 0.26 dex, which is 
expected simply from the metallicity variation among the four 
galaxies and within each galaxy (e.g., Willner & Nelson-Patel 
2002). 

Given the relation between the [Ne II] + [Ne III] luminos- 
ity and the ionizing photon rate, it is straightforward to trans- 
late L[Ne ii]+[Ne in] into a SFR, using one of the many existing 
calibrations between SFR and N- ion . Adopting the calibration 
of Kennicutt (1998), which assumes solar abundances and a 
Salpeter (1955) initial mass function with a lower mass limit 
of 0. 1 Mq and an upper mass limit of 100 M Q , and choosing a 
fiducial value of fi on = 0.6, 



SFR (Mq yr _1 ) = 4.34x 10 



-41 



L[Ne II]+[Ne III] (ergS S ) 



/++1-67/+2 



(13) 



£[Ne ii]+[Ne in] (ergs s 1 ) = 3.16 x 10 42 

The Ne luminosity can be easily related to the Lyman con- 
tinuum luminosity by recognizing that, in equilibrium, the ion- 
izing rate of hydrogen is equal to the recombination rate: 

Nioo (photons s _I ) = n e ii v a^V = «b EM A, (11) 

where as = 2.6 x 10~ 13 cm 3 s" 1 is the Case B recombination 
rate and V is the volume of the emitting region. Combining 
Equations 10 and 11, 

L[Ne ii]+[Ne in] (ergs s" 1 ) = 4.15 x 10" 13 

x CA + 1.67/ +2 ) / ion N ion (photons s" 1 ). (12) 

Here, f mn is the number fraction of ionizing photons that are 
actually absorbed by the gas (as opposed to being absorbed by 
dust or freely escaped from the nebula). Hirashita et al. (2003) 
estimate that /;„„ = 0.6 ± 0.2 in starburst galaxies. 

To check the validity of Equation 12, we present in Figure 4 
a compilation of [Ne II], [Ne III], and Bra measurements for a 
large collection of H II regions observed with the ISO SWS. The 
data for the H II regions in the Galaxy and the Small and Large 

1 Giveon et al. (2002) also give data for some M33 H II regions, but the data 



4. DISCUSSION AND SUMMARY 

We have investigated the utility of the mid-IR fine-structure 
lines of Ne as a SFR indicator. Being dominant coolants in 
H II regions, the [Ne II] and [Ne III] lines are among the most 
prominent and easy-to-measure features in the mid-IR spectra 
of virtually all extragalactic sources, especially with the advent 
of Spitzer. Unlike commonly used SFR indicators such as the 
synchrotron radio continuum or even the IR luminosity, which 
relate to massive stars in a complex manner, the formation of 
the [Ne II] and [Ne III] lines is particularly simple. Their in- 
tensity scales with the emission measure of an ionized nebula, 
and hence directly with the Lyman continuum luminosity. This 
property makes them an extremely attractive tool to estimate 
the SFR in extragalactic systems. The obvious advantage of 
[Ne II] and [Ne III] over optical lines such as Ha and [O II] is 
their greatly reduced sensitivity to dust extinction. 

Using a sample of star-forming galaxies with available mid- 
IR spectra from the literature, we show that the [Ne II] luminos- 
ity correlates linearly with the total IR luminosity, over 5 orders 
of magnitude in luminosity and ranging from low-metallicity 
dwarf galaxies to high-metallicity, ultraluminous IR galaxies. 
The observed scatter in this correlation is only ~ 0.5 dex, which 
is remarkably tight, considering several likely sources of vari- 
ation. First, the neon abundances span at least a factor of 5 
(e.g., Verma et al. 2003). Second, aperture effects must con- 
tribute to at least part of the scatter. Both the SWS aperture and 
the effective aperture of the IRS (either in mapping or staring 
mode) are significantly smaller than the IRAS beam. Although 
one usually expects the IR emission to be centrally peaked in 

of Willner & Nelson-Patel (2002) is more extensive. 
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starburst systems, given the diversity of our sample it is quite 
likely that the IR emission arises from a larger area than is sam- 
pled by the [Ne II] measurements. The data set from Thornley 
et al. (2000) attempted to correct the IR emission for aper- 
ture effects, but in general this was not available for the other 
data sets. To estimate the likely impact of aperture effects, we 
sorted the sample by optical angular diameter, but find no obvi- 
ous dependence of angular size scatter with scatter around the 
correlation. Third, the fractional ionization of Ne is variable. 
However, inclusion of both the [Ne II] and [Ne III] luminosi- 
ties improves the empirical correlation in that the scatter and 
turnover of the correlation at the low-luminosity end decrease 
with the combined luminosities. Lastly, while the extinction at 
mid-IR wavelengths is much less than in the optical, in some 
cases it can still be non-negligible (e.g., Genzel et al. 1998). 

A quantitative explanation of the L [Ne n]+[Ne ni] vs. L 1R corre- 
lation is nontrivial, since the IR emission arises not only from 
dust grain heating by massive stars but also by nonionizing 
stars, both young and old (e.g., Hirashita et al. 2003). How- 
ever, to a first approximation this correlation indicates that, to 
the extent that the IR luminosity traces the SFR (e.g., Kenni- 
cutt 1998; Bell 2003; Hirashita et al. 2003), so, too, does the 
[Ne II] + [Ne III] luminosity, as expected from the basic theo- 
retical considerations discussed in § 3. 

To our knowledge, this is the first systematic analysis of the 
[Ne II] and [Ne III] lines for a large sample of star-forming 
galaxies with the expressed intent of using them to derive SFRs. 
In their ISO survey of ultraluminous IR galaxies, Genzel et al. 
(1998) had explicitly recognized that the [Ne II] luminosity can 
be used as a surrogate for the Lyman continuum luminosity. 
From a subset of 12 starburst galaxies, they empirically deter- 
mined that the average ratio of Lyman continuum luminosity 
to [Ne II] luminosity is 64. This value is smaller than that de- 
duced from our calculations in §3: according to the starburst 
synthesis models presented in Hirashita et al. (2003), the av- 
erage Lyman continuum photon has an energy of <~ 19.4 eV, 
and if /+ = 0.75, f+2 = 0, and f on = 0.6, then Equation 12 
implies (L Lyc /L[ Ne irj) ~ 166. In a recent Spitzer IRS study 
of low-redshift quasars and starburst-dominated ultraluminous 
IR galaxies, Schweitzer et al. (2006) presented a correlation 
between [Ne II] luminosity and 60 /im luminosity, showing 
that the two quantities are significantly correlated. Although 



Schweitzer et al. argue that the [Ne II] emission in their sam- 
ple largely originates from star-forming regions rather than the 
accretion-powered narrow-line region, this is difficult to prove 
in the absence of a detailed analysis of the narrow-line spectra 
of their sample. In any event, our analysis spans a larger dy- 
namic range in luminosity and galaxy environment, and it con- 
clusively demonstrates that in star-forming galaxies the com- 
bined [Ne II] + [Ne III] luminosity tracks the IR emission. 

We propose that the relationship between the [Ne II] + 
[Ne III] luminosity and the SFR (Eq. 13) be added to the arsenal 
of SFR indicators presently used in the literature. If possible, 
the [Ne II] and [Ne III] luminosities employed in Equation 13 
should be corrected for extinction. We note, however, that our 
calibration itself, which is derived from theoretical considera- 
tions, does not depend on extinction correction. Our calibration 
is based on (1) the relation between the luminosity of the neon 
lines and the ionizing photon rate (Eq. 12) derived in this pa- 
per and (2) a previously determined and widely accepted con- 
version between the ionizing photon rate and the SFR (Kenni- 
cutt 1998). While the neon-based SFR calibration explicitly de- 
pends on the fractional abundance of singly and doubly ionized 
Ne (as expressed by the factors /+ and f+2 in Eq. 13), which 
can only be ascertained from detailed modeling of individual 
sources, if most of the Ne is ionized, the distribution between 
the [Ne II] and [Ne III] states introduces a variation in the cali- 
bration of less than a factor of 2. 

As an independent check of the reliability of our new method, 
we have compared the SFR obtained from the [Ne II] + [Ne III] 
luminosity with the SFR derived more traditionally from the IR 
luminosity for the sample of 57 galaxies plotted in Figure 1. 
Using the IR calibration given in Kennicutt (1998), the ratio of 
the IR-based SFRs to the Ne-based SFRs has an average value 
of 1.9 ± 3.7 and a median value of 0.7. Within the significant 
scatter, there is no evidence for a large systematic bias in the 
Ne-based method. 

This work was supported by the Carnegie Institution of 
Washington and by NASA grants from the Space Telescope 
Science Institute (operated by AURA, Inc., under NASA con- 
tract NAS5-26555). We thank Alex Dalgarno for discussions 
on collision rates. We are grateful to an anonymous referee for 
a thoughful and constructive review of the paper. 
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log (J-ib / L@) 

FIG. 1 . — Empirical correlation between ionized neon emission and the total IR (8-1000 /im) luminosity, derived using 12, 25, 60, and 100 pm flux densities from 
the Infrared Astronomical Satellite (IRAS) following the prescription of Sanders & Mirabel (1996). The open points denote [Ne II] 12.8 pm measurements, whereas 
the solid points represent the sum of [Ne II] and [Ne III] 15.6 fim, with the linear regression fits shown as dashed and solid lines, respectively. The measurements 
were taken from Genzel et al. (1998; excluding the active galactic nuclei), Thornley et al. (2000), Verma et al. (2003), O'Halloran et al. (2006), and Wu et al. 
(2006). We have used the distances adopted by these authors. Data taken with ISO SWS are plotted as circles; data taken with Spitzer IRS are plotted as squares; 
five sources have upper limits on Lir. If both ISO and Spitzer data are available, preference is given to the latter. 
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log (u / Hz) 

FIG. 2. — Photoionization cross sections for Ne° (solid, black), Ne + (dashed, red), and Ne +2 (dot-dashed, blue). The cross sections are computed using Equation 
2.31 and Table 2.7 of Osterbrock (1989). Shown for comparison are curves for blackbody emission at 4 X 10 4 and 6 X 10 4 K (dotted), which have been scaled 
arbitrarily. 
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FIG. 3. — Fractional abundance of Ne° (solid, black), Ne + (dashed, red), and Ne +2 (dot-dashed, blue) as a function of (a) radiation temperature for two fixed 
values of the ratio of the dilution factor Y for the blackbody radiation and the density n e (T/n e = 10~' 9 , right set of curves; 10~ 18 cm 3 , left set of curves) and (b) 
the ratio r/n c for a fixed radiation temperature of T r = 6 X 10 4 K. The calculations utilize the photoionization cross sections from Figure 2, radiative recombination 
rates from Gould (1978), and dielectronic recombination rates from Nussbaumer & Storey (1987). The recombination rates assume a fiducial temperature of 10 4 K 
and that the population of each species is in the ground state. The photoionization rate of [Ne II] includes transitions to the metastable states of [Ne III] ( l D and '$). 
We have ignored charge-exchange transitions involving the exchange of ionization state between neutral H and the Ne ions, since the density of neutral H should be 
negligible in ionized gas. 
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FIG. 4. — Empirical correlation between ionized neon emission and the Bra emission for H II regions (Giveon et al. 2002; Willner & Nelson-Patel 2002) in the 
Galaxy (circles), the Large Magellanic Cloud (squares), the Small Magellanic Cloud (triangles), and M33 (pentagons). The Br/3 fluxes for M33 were translated to 
Bra assuming a Bracket decrement of 1.6. The solid line shows the relation given in Equation 12, expressed in terms of Bra using the effective recombination rate 
for Bra for Case B recombination (rc c = 10 2 cnr 3 , T e = 10 4 K; Osterbrock 1989) and assuming /+ = 0.75 and f\ on = 0.6. 



